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Report  Title 

Micromechanics  Based  Representative  Volume  Element  Modeling  of  Heterogeneous  Cement  Paste 

ABSTRACT 

The  current  work  focuses  on  evaluation  of  the  effective  elastic  properties  of  cementitious  materials  through  a  voxel 
based  FEA  approach.  Voxels  are  generated  for  a  heterogeneous  cementitious  material  (Type-I  cement)  consisting  of 
typical  volume  fractions  of  various  constituent  phases  from  digital  microstructures.  The  microstructure  is  modeled  as 
a  micro-scale  representative  volume  element  (RVE)  in  ABAQUS  to  generate  cubes  several  tens  of  microns  in 
dimension  and  subjected  to  various  prescribed  deformation  modes  to  generate  the  effective  elastic  tensor  of  the 
material.  The  RVE-calculated  elastic  properties  such  as  moduli  and  Poisson’s  ratio  are  validated  through  an 
asymptotic  expansion  homogenization  (AEH)  and  compared  with  rule  of  mixtures.  Both  Periodic  (PBC)  and 
Kinematic  boundary  conditions  (KBC)  are  investigated  to  detennine  if  the  elastic  properties  are  invariant  due  to 
boundary  conditions.  In  addition  the  method  of  “Windowing”  was  used  to  assess  the  randomness  of  the  constituents 
and  to  validate  how  the  isotropic  elastic  properties  were  detennined.  The  average  elastic  properties  obtained  from  the 
displacement  based  FEA  of  various  locally  anisotropic  micro-size  cubes  extracted  from  an  RVE  of  size  100x100x100 
microns  showed  that  the  overall  RVE  response  was  fully  isotropic.  The  effects  of  domain  size,  degree  of  hydration, 
kinematic  and  periodic  boundary  conditions,  domain  sampling  techniques,  local  anisotropy,  particle  size  distribution 
(PSD),  and  random  microstructure  on  elastic  properties  are  studied. 
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Abstract 


□  Effective  elastic  properties  of  cementitious  materials  are  evaluated 
through  a  voxel  based  FEA  approach. 


Introduction 


□  A  methodology  has  been  developed  for  computing  the  elastic  properties 
of  heterogeneous  C-S-H  (calcium  oxide-  silicate  oxide-  hydroxide)  based 
multi-phase  cementitious  materials. 

□  The  primary  focus  is  to  predict  homogenized  properties  at  macro-levels 
using  micro  mechanics  based  models. 

□  Focus  is  on  the  determination  of  elastic  properties  for  hydrated  cement 
paste  from  un-hydrated  constituents  when  small  strain  quasi-static  loading 
conditions  are  applied  to  micro-scale. 
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Methodolo 


□  The  representative  volume  element  (RVE)  is  the  smallest  volume  of 
material  that  captures  global  characteristics  of  the  material  and  shows  the 
same  overall  material  properties  irrespective  of  the  boundary  conditions 
applied. 

□  Software  package  CEMHYD3D  V.3  (NIST),  simulates  the  hydration 
process  and  formation  of  the  digitally  generated  micro -structure  for  a 
typical  Type-I  general  purpose  cement. 

□  Initial  3D  micro  structure  is  created  based  on  measured  geometrical 
particle  size  distribution  (PSD)  as  well  as  volume  fractions  and  surface- 
area  fractions  of  the  constituent  phases  for  cement  powder,  extracted  from 
2D  composite  images  of  cement  at  various  degrees  of  hydration  (DOH). 

□  The  RVE-calculated  elastic  properties  such  as  moduli  and  Poisson’s  ratio 
are  validated  through  an  asymptotic  expansion  homogenization  (AEH) 
and  compared  with  rule  of  mixtures. 

□  Windowing  is  employed  to  investigate  how  anisotropy  due  to  local 
micro  structure  leads  to  overall  isotropic  behavior  of  the  agglomerate. 
Windows  are  analogous  to  physical  core  samples  prepared  by  extraction 
from  a  hydrated  bulk  specimen. 


genpartnew.c 


Input: 

•discretized  PSD  for  cement  clinker, 
gypsum,  slag,  calcium  carbonate,  inert 
filler 

•number  of  particles  to  place 
•proportion  of  “cement”  particles  that 
should  be  calcium  sulfate  (gypsum, 
hemihydrate,  and/or  anhydrite), 

•whether  particles  should  be  flocculated, 
random  or  dispersed  during  placement, 
•whether  aggregate  should  be  included 
(for  studying  ITZ) 


Output:  initial  3D 
microstructure  (periodic  B.C) 

•  microstructure  (Phase  ID)  voxel 
(cement  particle)  coordinates 

•  individual  particle  IDs  (for  assessing 
setting  of  cement  during  hydration) 


distrib3d.c 


Input: 

•  initial  3D  microstructure  (Phase  ID) 

•  measured  2D  correlation  files  for 
major  phase  combinations  (silicates, 
C3S.C3A,  C4AF) 


Output:  phases  assigned  to 
voxels 

•  Updated  microstructure  (Phase  ID) 
•four  major  cement  clinker  phases  (C3S, 
C2S,  C3A,  and  C4AF)  distributed 
amongst  voxels  in  3D  microstructure 
•user-requested  volume  and  surface 
area  fractions  for  each  phase  achieved 


stat3d.c 


Input:  3-Dimage 


Output:  Reports 

•phase  volumes 
•volume  fractions 
•pore-exposed  surface  area 
fractions 
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C3S-1,  C2S-2 ,  C3A-3,  C4AF-4, 
gypsum-5,  hemihydrate-6,  anhydrite- 
7,  aggregate-28,  and  CaC03-26 
(separated  by  spaces).  Typical 
values  following  phase  distribution 
using  the  distrib3d.c  program  would 
be,  2,  3, 4,  5,  6,  7,  26,  and  28. 


disrealnew.c 


Input: 

•3-D  microstructure  (Updated 
Phase  ID)  generated  by 
distrib3d.c 

•  particle  IDs  (for  assessing 
setting  of  cement  during 
hydration) 

•Phase  volume  fractions 
•Phase  area  fractions 


Output:  cement 
hydration 

•Uses  cellular  automata 
techniques 

•Preserves  correct  hydration 
volume  stoichiometry 


3D  Mesh  generation 
(brick  elements) 


Generation  of  particle  size  distributions,  voxel  coordinates  and  material  type 


Phase 

Phase 

ID 

Volume 

Fraction0/*) 

E 

(GPa) 

Poisson’s 

ratio 

V 

1 

Water  Porosity 

0 

19.8251 

0.001 

0.499924 

2 

Tricalcium  Silicate  (C3S) 

1 

7.1625 

117.6 

0.314 

3 

Dicalcium  Silicate  (C2S) 

2 

2.628 

117.6 

0.314 

4 

Tricalcium  Aluminate  (C3A) 

3 

1.7376 

117.6 

0.314 

5 

Tetracalcium  Aluminoferrite  (C4AF) 

4 

1.1012 

117.6 

0.314 

6 

Dihydrate  (Gypsum)  (('.S’.  H? ) 

5 

0.0022 

45.7 

0.33 

7 

Hemihydrate  ( CS.H1/2  ) 

6 

0.0001 

62.9 

0.3 

8 

Anhydrite  (CS) 

7 

0.0005 

80 

0.275 

9 

Calcium  Hydroxide  (CH) 

13 

14.425 

42.3 

0.324 

10 

Calcium  Silicate  Hydrate  Gel  (CSH) 

14 

37.4425 

22.4 

0.25 

11 

Hydrogarnet  (C3AH6) 

15 

4.2538 

22.4 

0.25 

12 

Ettringite  (QASAC?) 

16 

6.034 

22.4 

0.25 

13 

Iron-rich  Stable  Ettringite 
(ETTRC4AF) 

17 

1.807 

22.4 

0.25 

14 

Monosulfate  AFM(C4ASH12  ) 

18 

2.4623 

22.4 

0.25 

15 

Iron  Hydroxide  (FH3) 

19 

0.3193 

22.4 

0.25 

16 

Gypsum  Formed  from  Hemihydrate 
and  Anhydrite  (GYPSUMS) 

25 

0.003 

45.7 

0.33 

17 

ABSGYPS 

29 

0.2996 

45.7 

0.33 

18 

Empty  Porosity 

45 

0.4963 

0 

0 

Material  properties  and  volume  fractions  of  constituent  phases 
for  a  representative  Type-I  cement 

Linear  elastic 
Nonlinear  Quasi-static 
Nonlinear  Dynamic  Explicit 
(wave  propagation) 


10000 


1000 


100x100x100 

50x50x50 

-20x20x20 

10x10x10 


4  6  8  10  12  14  16  18 

Particle  Radius  (micron) 


Scaling  is  required  to  maintain  consistent  area  and 
volume  fractions  of  various  constituent  phases 


Cement  Hydration 


Hydrated  microstructure 
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Hydration  Products 
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Typical  volume  fractions  of  major  constituents  vs  Curing 
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Development  of  the  Young’s  Modulus 
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Volume  fractions  of  major  phases  for  IK,  8K,  125K  and  1M  PMDs 
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Constraint  Equations  for  PBC 
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□  Subject  to  PBC,  exact  estimate  of  the  effective  homogeneous  elastic  properties  can  be  obtained 
for  linear  elastic  inhomogeneous  microstructures  that  exhibit  perfectly -periodic  homogeneity  by 

77111 

solving  lor Xk  in: 

□  Vector  Vj  signifies  the  coordinates  of  the  microstructure  RVE,  and  Dijkl  is  the  elastic  stiffness 

^  1 

tensor  at  a  point  y  in  the  material.  The  homogenized  linear  elastic  stiffness  tensor,  D  is  given 


Microstructure  Based  Homogenization 


Effect  of  DOH  on  material  bulk  properties  for  1M-RVE  for  KBC  and  PBC 
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Effect  of  domain  size  on  material  bulk  properties  for  DOH=0.8  for  KBC  and  PBC 
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Volume  fractions  of  major  phases  for  (a)  IK  element  and  (b) 

8K  element  windows 


Elastic  moduli  (uniaxial)  with 
increasing  window  size 


Elastic  shear  modulus  Gij  with 
increasing  window  size 
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□ 


□ 


Rule  of  Mixtures  Based  Homogenization 


A  rule  of  mixtures  approach  independent  of  the  microstructure  of  the  material  is  used 
to  compute  the  effective  bulk  properties  of  the  cementitious  material. 

The  theoretical  extreme  upper  and  lower  bounds  on  effective  material  properties  of 
multi-phase  materials  are  the  Voigt  (1928)  and  Reuss  (1929)  bounds. 
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Hashin(1962)  presented  the  composite  (or  coated)  spheres  model  for  determining  the 
effective  material  properties  for  multi-phase  materials,  based  on  the  dilute  suspension 
model. 
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□  For  Hashin  and  Voigt  estimates,  the  bulk  modulus  (K)  is  found  to  be  lower  compared 
to  the  values  computed  based  on  the  microstructure  (KBC,  PBC  and  AEH).  Both  the 
Young’s  Modulus  (E)  and  shear  modulus  (G)  are  determined  to  be  higher  than  those 
estimated  by  microstructure  based  homogenization. 


Conclusion 


□  A  comparison  between  the  two  domain  sampling  methods  shows  that  windowing  produces 
effective  material  properties  with  a  larger  variation  than  the  PMD  due  to  a  higher  variation 
in  local  phase  volume  fractions. 

□  Macroscopic  properties  obtained  for  various  DOH  and  domain  sizes,  determined  by 
applying  Kinematic  Boundary  Conditions  (KBC),  Periodic  Boundary  Conditions  (PBC), 
AEH  and  rule  of  mixtures  based  homogenization  are  found  to  be  comparable. 

□  It  is  shown  that  even  though  cement  is  a  heterogeneous  anisotropic  material  at  the  micro¬ 
level,  the  bulk  properties  are  effectively  isotropic. 
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